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Abstract. - We investigate the transition between the Cassie-Baxter and Wenzel states of a 
slowly evaporating, micron-scale drop on a superhydrophobic surface. In two dimensions analytical 
results show that there are two collapse mechanisms. For long posts the drop collapses when it is 
able to overcome the free energy barrier presented by the hydrophobic posts. For short posts, as 
the drop loses volume, its curvature increases allowing it to touch the surface below the posts. We 
emphasise the importance of the contact line retreating across the surface as the drop becomes 
smaller: this often preempts the collapse. In a quasi-three dimensional simulation we find similar 
behaviour, with the additional feature that the drop can de-pin from all but the peripheral posts, 
so that its base resembles an inverted bowl. 
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Introduction. — It is well-known that the hydropho- 
bic nature of a surface is amplified by its roughness [1,2]. 
This can happen in two different ways. When the liquid 
drop occupies the spaces between the surface projections, 
and is everywhere in contact with the surface, it is said to 
Be in the collapsed or Wenzel state [3] . The contact angle 
is 



cos 6 



W 



(1) 



where r is the ratio between the real surface area and its 
projection onto the horizontal plane and 9 e is the equilib- 
rium contact angle of the flat surface. On the other hand, 
if penetration does not occur and the drop remains bal- 
anced on the surface projections with air beneath it, it is 
in the suspended or Cassie-Baxter state [4] with contact 
angle 



cosOcb = $>cos9 e — (1 - $) 



(2) 



with $ the solid fraction of the surface. Both states are 
(local) minimum of the free energy, but there is often a fi- 
nite energy barrier opposing the transition between them. 
The magnitude of the energy barrier has been shown to 
depend on both the size of the drop and the roughness of 
the surface [5,6]. 

The main aim of this paper is to explore the mecha- 
nisms by which the drop spontaneously collapses [7, 8] . 
We consider micron-scale drops, sufficiently large that we 



can ignore thermal fluctuations but smaller than the cap- 
illary length so that gravity is not important. We focus on 
the limit where the evaporation timcscale is much longer 
than the timcscale for drop equilibration so that the drop 
is always in thermodynamic equilibrium. This is normally 
the physically relevant situation for experiments on micron 
scale drops. The question of how and when collapse occurs 
is important because, even though both states show high 
values of the contact angle, many of their other physical 
properties, for example, contact angle hysteresis are very 
different [9]. 

We first consider a drop on a two dimensional, superhy- 
drophobic surface and present analytic results for how it 
collapses as its volume is decreased. We argue that there 
are two mechanisms for collapse. For short posts, as the 
curvature of the drop increases, it touches the surface be- 
low the posts, thus breaching the free energy barrier. For 
longer posts the free energy barrier is removed when the 
surface free energy gained by the drop as it collapses wins 
over the surface free energy lost by increased contact with 
the hydrophobic posts. However, importantly, the collapse 
transition is usually preempted by the contact line of the 
drop retreating across the surface. Therefore collapse for 
drops on long posts will normally occur only when the 
drop covers a very small number of the posts. 

In three dimensions anaytical calculations are not feasi- 
ble so we use numerical simulations to follow the behaviour 
of the shrinking drop. A new feature is that the base of the 



H. Kusumaatmaja et al. 




Fig. 1: Schematic diagram of a suspended drop. 



drop tends to form a bowl-shape, where the lines of con- 
tact depin and move down all but the outermost posts. 
We further argue that the tendency for the contact line 
to prefer to retreat across the surface than to collapse is 
even more pronounced in three dimensions than in two. A 
conclusion summarises our results and compares them to 
experiments. 

Drop collapse in two dimensions: analytical re- 
sults. — Consider a two dimensional drop suspended on 
a regular array of hydrophobic posts as shown in Fig. 1. 
The posts have width a, spacing b and height I, and the 
substrate material has an intrinsic contact angle 9 e > 90°. 
The drop forms a circular cap with a contact angle 9, cross- 
sectional area S, radius of curvature R and base length 2r. 

We consider a drop with contact line that is pinned at 
the outer edges of two posts and we first assume that the 
contact line pinning persists as the drop spontaneously col- 
lapses. Considering the motion of a retreating contact line 
across a superhydrophobic surface (as, say, the volume of 
the drop is slowly decreased) shows that the line is pinned 
for 180° > 8 > 6 e [9]. After we have explained the possi- 
ble collapse mechanisms we shall return to a discussion of 
when they are preempted by dcpinning. 

Given pinning, the base radius r is fixed and takes dis- 
crete values 

r = (m+ 1/2) a + m 6 (3) 

where 2m + 1 = 1, 2, 3, ... is the number of posts beneath 
the drop. The cross-sectional area of the drop, which is 
constant, can be written 



S = 



— sin v cos ( 
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The last term in Eq. (4) is due to the curved interface 
underneath the drop and </> = 9 p — 90° where 9 p is the 
angle this interface makes with the sides of the posts. 

Our aim is to investigate when and how the collapse 
transition occurs. We do this by considering the behaviour 
of the drop free energy as a function of h, the distance 
it penetrates into the substrate (see Fig. 1). The non- 
constant contributions to the drop free energy F come 



from three terms. The first two correspond to the liquid- 
gas interfacial free energy above and beneath the surface 
and the third term is the free energy required by the liquid 
drop to wet the posts to a depth h 



f = Fh = 



2r0 
sin 



2 m b <p 
sin d> 



Amh cos 9 P 



(5) 



where 7 is the liquid-gas interfacial tension. 

We now consider the variation of the free energy with 
h. The drop will start to penetrate the posts if ^ < at 
h = 0, or equivalently % > 0, since 4| < 0. Using the 
constraint of constant area to eliminate dh gives 

2r (sin 6 -9 cos 9) , . . 2r n . Jn 

df = — 5 (sm 9 + — cos 9 e ) d9 + 

sin 9 b 



2mb (sin <p — (f> cos < 



(sine/) + cos 9 e ) d(j> ■ (6) 



Consider first d<j> = 0. Since 2 r (sin 9 — 9 cos 9) / sin 3 9 > 
0, the condition for the drop to start collapsing is 



2r 



■ cos 9 e > 



h=0 



(7) 



The corresponding critical drop radius of curvature and 
contact angle are [10] 



Rr = —. 



sin 9 r 



f cos 9 e . 



(8) 
(9) 



9 gets smaller and sin 9 gets larger as the drop penetrates 
the posts. As a result, once Eq. (7) is satisfied it will 
always be satisfied and once the drop has started to move 
it collapses fully, into the Wenzel state. 

The drop will be in equilibrium at h = on the thresh- 
old of the collapse transition. Therefore we may combine 
Eq. (9) and the Laplace pressure condition to show that 
9 p = 9 e , or <f> = 6 e — 90° as expected from the Gibbs' 
criterion [11]. Hence, from Eq. (6) the free energy is at 
an extrcmum with respect to changes in Calculating 
the second derivative confirms that this is a minimum and 
hence that the assumption d<f> = is appropriate. 

Typical plots of the free energy of a drop against h, the 
distance it penetrates into the substrate are shown in Fig. 
2, where for simplicity we have neglected the corrections 
due to the curvature of the interfaces in the grooves. In 
Fig. 2(a), where we have used m = 3, b/a = 1.5, 9 e = 
95°, and 9\ h=0 = 111° < 9 C = 111.6° the free energy 
is a smoothly decreasing function of h and the drop will 
collapse. In Fig. 2(b) on the other hand, for 9\h=o = 
112° > 9 C = 111.6°, there is a free energy barrier and 
therefore no collapse. The peak of the free energy barrier 
occurs at 9 = 9 C and the magnitude of the barrier is 
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Fig. 2: Normalised drop free energy against penetration depth 
when (a) the collapse transition occurs and (b) there is a free 
energy barrier between Cassie-Baxter and Wenzel states. Fo 
is the drop free energy in the Cassie-Baxter state, 7 is the 
liquid-gas surface tension and S is the drop area, m = 3, 
b/a = 1.5, 9 e = 95° and 9\ h =o = 111° and 112° for (a) and (b) 
respectively. 



We have argued that, for R < R c , there is no free en- 
ergy barrier to drops penetrating hydrophobic posts. The 
critical radius depends on the post width a, the post sep- 
aration &, the base radius r, and the equilibrium contact 
angle 9 e . It does not, however, depend on the post height 
I. There is, however, another route to drop collapse [8], 
which will pre-empt this mechanism for shallow posts. 

Prior to collapse the liquid drop has not penetrated the 
posts, the system is in mechanical equilibrium, and the 
Laplace pressure is the same everywhere. Thus the liquid- 
gas interface between the posts bows out with a radius of 
curvature equal to that of the circular cap R. The centre 
of the curved interface reaches a distance d into the posts: 



8R 



(11) 



for small (j>. As R gets smaller, d increases. When d = I the 
liquid-gas interface touches the base surface initiating the 
transition between the Cassie-Baxter and Wenzel states. 
At this point there is a considerable free energy release 
because the drop is replacing two interfaces (liquid-gas 
and gas-solid) with a single liquid-solid interface. Conse- 
quently this transition is irreversible and for the opposite 
transition to occur (Wenzel to Cassie-Baxter) an external 
force is needed to overcome the free energy barrier. For 
this mechanism to be possible it is apparent from simple 
geometry that d < 6/2. 

Regions of parameter space where there is (i) collapse 
due to the contact line sliding down the posts, (ii) collapse 
due to the centre of the interface touching the base surface, 
(iii) no collapse are distiguished in Fig. 3. The crossover 
between regions (i) and (ii) occurs when 

cos6» e < — -. (12) 
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Fig. 3: The crossover between the two different drop collapse 
mechanisms in two dimensions. 



It is interesting to note that the crossover point between 
the two regimes (Eq. 12) will slide to larger I as the posts 
are made more hydrophobic. 

We now revisit the assumption that the contact line is 
pinned at the outer edges of the posts. The (theoretical) 
advancing contact angle is 180° [9] and therefore the line 
will not move outwards. The receding angle in the quasi- 
static limit is 9 e [9] and therefore it will not jump inwards 
if 9 C > 9 e or, equivalently, sin 2 9 C < 1 — cos 2 6> e . Using 
Eqs. (3) and (9) this is equivalent to 



cos 2 9 e < 



. m + 1/2 
4 I — — -!— + m 

b/a 



+ 1 



(13) 



Fig 4(c) shows the maximum value of 9 e at which col- 
lapse will occur for different m and b/a. For small b/a col- 
lapse is strongly suppressed and only occurs for tiny drops 
on slightly hydrophobic surfaces. Even for b/a >> 1 the 
tendency to depin is strong and the collapse occurs for 
small value of m unless 9 e is close to 90°. By setting r = b 
in Eq. (7) , we conclude that the transition will never occur 
spontaneously for 9 e > 120°. 

In Fig. 4(a) and (b), we plot the critical radius of cur- 
vature at which a transition occurs R c /a and the corre- 
sponding drop area (which we present as \J S/ira 2 ) as a 
function of 9 e , m and b/a. As expected the critical radius 
of curvature does not depend on m; it is a function of b/a 
and 9 e only. The critical base area of the drop does, how- 
ever, depend on m and is smaller for larger values of m 
and a/b. The curves for increasing m terminate at decreas- 
ing values of 9 e corresponding to the contact line receding 
inwards before the drop is able to penetrate the posts. 

Simulations of drop collapse. We now describe 
the details of a numerical model which will allow us to 
explore the collapse transition in both two and three di- 
mensions. We describe the equilibrium properties of the 
drop by a continuum free energy [12] 



(Mn) + ~{d a nf)dV + 



Mn s )dS. (14) 
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Fig. 4: (a) Critical drop radius of curvature and (b) area (presented as */ S/jva 2 ) for collapse as a function of the equilibrium 
contact angle 8 e for different substrate geometries and number of posts beneath the drop, (c) The maximum value of 6 e for 
which a spontaneous drop collapse can occur for different b/a and m. 
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Fig. 5: Critical drop area (presented as \J S/iva 2 ) as a function 
of contact angle: comparison between the two-dimensional the- 
ory and simulations for m = 3, a = 8, b = 12. The circles and 
crosses represent drops which collapse and remain suspended 
respectively. The solid line is the theoretical prediction. 



ipb(n>) is a bulk free energy term which we take to be [12] 
Mn) = Vc(y n + 1)>' - 2f„ + 3 - 2(3t w ) , (15) 

where v n = (n — n c )/n c , t w = (T c -T)/T c and n, n c , 
T, T c and p c are the local density, critical density, local 
temperature, critical temperature and critical pressure of 
the fluid respectively. This choice of free energy leads to 
two coexisting bulk phases of density n c {H /3t w ) , which 
represent the liquid drop and surrounding gas respectively. 
Varying (3 has the effects of varying the densities, surface 
tension, and interface width; we typically choose (3 = 0.1. 

The second term in Eq. (14) models the free energy 
associated with any interfaces in the system, k is related 
to the liquid-gas surface tension and interface width via 
a ig = (4V2^(/3r tu ) 3 / 2 n c )/3 and £ = {nnl/A(3T wPc ) 1 / 2 
[12]. We use k = 0.0018, p c = 1/8, t w = 0.3, and n c = 3.5. 

The last term in Eq. (14) describes the interactions be- 
tween the fluid and the solid surface. Following Calm [13] 



the surface energy density is taken to be ip s (n) = —Xn s , 
where n s is the value of the fluid density at the surface. 
The strength of interaction, and hence the local equilib- 
rium contact angle, is parameterised by the variable A. 
Minimising the free energy (4) leads to a boundary condi- 
tion at the surface, d±n = —X/k, and a relation between 
A and the equilibrium contact angle e [12] 



A = 2(3r w ^2p c K sign(- - 9 e )]J cos -(1 - cos-) , (16) 

where a = cos" 1 (sin 2 9 e ) and the function sign returns 
the sign of its argument. Similar boundary conditions can 
be used for surfaces that are not flat: a way to treat the 
corners and ridges needed to model supcrhydrophobic sur- 
faces is described in [14]. 

The equations of motion of the drop are the continuity 
and the Navier-Stokes equations 

d t n + d a (nu a ) = , (17) 
d f (nu a ) + d/3(nu a up) = -dpPap + 
vdp [n(dpu a + d a up + S a pd 7 u-y)] (18) 

where u, P, and v are the local velocity, pressure tensor, 
and kinematic viscosity respectively. The thermodynamic 
properties of the drop appear in the equations of motion 
through the pressure tensor P which can be calculated 
from the free energy [12, 14] 

P&P = ipb ~ f {d a n) 2 - Knd^n)S al3 + n(d a n)(dpri), 

Pb = Pc(v n + 1) 2 (3^ 2 - 2v n + 1 - 2/3t w ). (19) 

When the drop is at rest d a P a p = and the free energy 
(14) is minimised. As we are considering the quasi-static 
problem when the drop is in equilibrium until the point 
of collapse details of its dynamics should not affect the 
results. However we choose to implement physical equa- 
tions of motion as this helps the drop to reach equilibrium 
quickly as its volume is decreased and for comparison to 
possible work on non-equilibrium collapse. 
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Fig. 6: (Color online) Equilibrium drop configurations, (a) 
9 e — 95° and varying volume, (b) Fixed volume and contact 
angles, 9 = 93° (black), 95° (purple), 97° (blue), 100° (green), 
and 110° (red). 



We use a lattice Boltzmann algorithm to solve Eqs. (17) 
and (18). No-slip boundary conditions on the velocity 
are imposed on the surfaces adjacent to and opposite the 
drop and periodic boundary conditions are used in the two 
perpendicular directions. Details of the lattice Boltzmann 
approach and of its application to drop dynamics are given 
in [12,14-17]. 

To implement evaporation we need to slowly decrease 
the drop volume. To do this we vary the liquid density by 
—0.1% every 2 x 10 5 time steps to ensure that the evapo- 
ration timescale is well separated from the drop equilibra- 
tion timcscalc. This in turn affects the drop volume as the 
system relaxes back to its coexisting equilibrium densities. 

Results for two dimensions are compared to the analytic 
solution in Fig. 5, where we have used a = 8, b = 12, and 
e = 95°. The critical drop area at which the collapse 
transition occurs is close to the theoretical value but criti- 
cal radii obtained from simulations arc typically too large 
by ~ 2 lattice spacings. This is because the liquid-gas 
interface is diffuse (~ 3 — 4 lattice spacings). We checked 
that, as expected, R c is independent of the post height 
and that h is the same everywhere underneath the drop. 

Drop collapse in three dimensions: numerical re- 
sults. — Analytic calculations in three dimensions are, 
in general, not possible for several reasons. Firstly, the 
drop shape is not a spherical cap but is influenced by the 
underlying topological patterning. Secondly, the shape of 
the liquid-gas interface spanning the posts is complicated. 
Thirdly, h, the distance the drop penetrates the substrate, 
is not neccesarily the same everywhere. Therefore we need 
to use the numerical approach presented in the last sec- 
tion to explore collapse. We consider a square array of 
posts of widths a = 3 and spacing 6 = 9. We present re- 
sults for both spherical drops and 'cylindrical' drops which 
demonstrate the relevant physics but are less demanding 
in computer time. 

A new feature in three dimensions is that for a spherical 
(or cylindrical) drop on a square array of posts the base of 
the drop can form a bowl-shape where the lines of contact 
with the top of all but the peripheral posts depin and 
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Fig. 7: (Color online) Evolution with time of a cylindrical drop 
on a square array of posts of width a = 3, spacing b — 9 and 
height I = 15. (a-c) Time evolution before collapse showing 
depinning of the receding contact line (note the scale change 
between (b) and (c)). (d— f) Motion of the collapsing drop: (d) 
cross sections in the plane bisecting the posts, (e) same times 
as (d), but in the plane bisecting the gap between the posts, (f) 
cross sections in the plane bisecting the gap, but with / = 45 
to enable the collapse to be followed to later times. 



move down the posts leaving the drop suspended by just 
its outer rim. This was seen in simulations for drops of 
both cylindrical and spherical symmetry, and has recently 
been reported experimentally [19]. The depinning occurs 
to reduce the distortion of the interface from spherical. 
Depinning is favoured for smaller drops and for contact 
angles close to 90°. 

This is apparent in Fig. 6(a) which shows the equi- 
librium profile of the drop as its volume is varied with 
B e = 95°. As expected the drop penetrates further into 
the posts as the radius is decreased (corresponding to in- 
creasing curvature). Fig. 6(b) shows cross sections of the 
final states of 5 drops with volumes V — 4.5 x 10 5 (in the 
units of lattice spacing) equal to within ~ 1% but varying 
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equilibrium contact angles in the range 9 e = 93° to 110°. 
Although the drop penetrates deeper into the posts as the 
intrinsic contact angle approaches 90° there is no collapse 
(these values would give a collapse transition in two di- 
mensions). Instead, for 9 e — 110°, the contact line depins 
and moves to cover 9 rather than 21 posts. Note that this 
jump also corresponds to a transition to the state when 
the drop is suspended on all the posts beneath it. not just 
those around its rim. 

To explore the depinning further, and to try to find a 
collapse transition, we turned to the geometry of a cylin- 
drical drop on a square array of posts. This preserves 
the physics whilst allowing us to exploit the quasi-two- 
dimensional geometry to run larger simulations. Results 
are shown in Fig. 7 for 9 e = 93°. Successive frames 
show how the drop profile evolves as its volume is quasi- 
statically decreased (note that they are drawn on different 
scales). Initially, the contact line is pinned at the edges 
of the posts and the drop penetrates further beneath the 
posts as the radius is decreased. However, as the drop con- 
tinues to decrease in size, the drop contact angle reaches 
the receding angle and the contact line depins. As it de- 
pins we observe that the penetration into the posts de- 
creases (because the drop is approximately spherical and 
the base area is reduced), thus moving the system away 
from the point where either a curvature or a free energy 
driven collapse is favourable. Eventually collapse is seen 
but only, for this example, when the drop spans just three 
posts. Note that for I — 45 the drop stops moving once 
it is fully inside the posts as its free energy becomes in- 
dependent of height: it forms a liquid bridge connecting 
several neighbouring posts. 

Indeed, we expect from the two dimensional calcula- 
tions that collapse is preempted by depinning for posts 
with b/a ~ 1 until the drops are very small. In three di- 
mensions depinning will be even more important because 
the receding contact angle is larger than 9 e , its value in two 
dimensions, because the distortion of the interface makes 
it more favourable for the drop to depin. From Fig. 7(a) 
and (b), we obtain 9 R ~ 120°. 

Summary. — To conclude, we have investigated the 
behaviour of an evaporating drop on a superhydrophobic 
surface. As the drop volume decreases quasi-statically it 
can move in three ways: (i) the drop attains its receding 
contact angle and the contact line moves inwards across 
the surface (ii) the free energy barrier to collapse van- 
ishes and the drop moves smoothly down the posts (iii) 
the drop touches the base of the surface patterning and 
immediately collapses. The depinning (i) is predominant 
and, unless the posts are widely spaced, or the surface 
is only very weakly hydrophobic, collapse occurs only for 
drops spanning a very small number of posts. 

This suggests strategies that could be used to suppress 
transitions to the Wenzel state. Long enough posts are 
needed to prevent curvature-driven collapse, i.e. / > b 2 /R, 
and the free energy barrier to the transition can be en- 



hanced by choosing 9 e as large as possible and using closely 
spaced posts, i.e. b < a. A mobile contact line will also 
help as this will relax any build up of curvature. 

Our results are in line with recent experiments [7,8,20]. 
In [8], for long posts, the contact line retreated as the 
drop shrank and collapsed only at the very end of evapora- 
tion. For short posts, a few depinning events were followed 
by collapse at a radius consistent with a curvature-driven 
mechanism, R c oc b 2 /l. It is not clear, however, whether 
the drop interface was suspended on all the posts or just 
those at the rim at the point of collapse: this detail is 
important in determining the constant of proportionality. 
In [19], the various drop configurations found here are also 
observed, including the depinning of the drop from all but 
the outer posts. In [20], for the somewhat different situ- 
ation of drops bounced onto a surface, the critical pres- 
sure for impalement varied linearly with post height for 
short posts, as expected for curvature-driven collapse, and 
showed a clear crossover to a length-independent regime 
for longer posts, consistent with a drop overcoming a free 
energy barrier. 

* * * 
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